arXiv:1504.00328vl [astro-ph.HE] 1 Apr 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 2 April 2015 (MN KTeX style file v2.2) 


Resonant oscillations of GeV - TeV neutrinos in internal shocks from 
gamma-ray burst jets inside the stars 


Nissim Fraija * 

Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, Circuito Exterior, 
C. U., A. Postal 70-264, 04510 Mexico D.F., Mexico 


2 April 2015 


ABSTRACT 

High-energy neutrinos generated in collimated jets inside the progenitors of gamma-ray bursts 
(GRBs) have been related with the events detected by IceCube. These neutrinos, produced 
by hadronic interactions of Fermi-accelerated protons with thermal photons and hadrons in 
internal shocks, are the only signature when jet has not broken out or failed. Taking into 
account that the photon field is thermalized at keV energies and the standard assumption 
that the magnetic field maintains a steady value throughout the shock region (with a width 
of 10^° — 10^^ cm in the observed frame), we study the effect of thermal and magnetized 
plasma generated in internal shocks on the neutrino oscillations. We calculate the neutrino 
effective potential generated by this plasma, the effects of the envelope of the star, and the 
vacuum on the path to Earth. By considering these three effects, the two (solar, atmospheric 
and accelerator parameters) and three neutrino mixing, we show that although GeV - TeV 
neutrinos can oscillate resonantly from one flavor to another, a nonsignificant deviation of the 
standard flavor ratio (1:1:1) could be expected on Earth. 
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1 INTRODUCTION 

Long gamma-ray bursts (IGRBs) have been associated to core col¬ 
lapse of massive stars leading to supemovae (CCSNe) of type Ib,c 
and II. Type Ic supernovae are believed to be He stars with radius 
Ri 10^^ cm, and type II and Ib are thought to have a radius 
of i?* ~ 3 X 10^^ cm. Depending on the luminosities and dura¬ 
tions, successful IGRBs have revealed a variety of GRB popula¬ 
tions: low-luminosity (11), ultra-long (ul) and high-luminosity (hi) 
GRBs (Meszaros & Waxman 2001; Liang et al. 2007; Gendre et al. 
2013). While llGRBs and ulGRBs have a typical duration of (~ 
10® - 10"* s), hlGRBs have a duration of tens to hundreds of sec¬ 
onds. Another important population associated with CCSNe, al¬ 
though unobservables in photons, are failed GRBs which could be 
much more frequent than successful ones, limited only by the ratio 
of type Ib/c and type II SNe to GRBs rates. This population has 
been characterized by having high-luminosities, mildly relativistic 
jets and durations from several to ten seconds (Huang et al. 2002; 
Meszaros & Waxman 2001; Soderberg & et al. 2006, 2010). 
Neutrinos are useful for studying the insides of stars, especially 
where photons cannot be observed either because jet fails or has 
not broken out yet, so in this case, they could be the only sig¬ 
nature that would display the dynamics of the star. High-energy 
(HE) neutrinos from this population of stars have been pointed out 
to contribute significantly to the extragalactic neutrino background 
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(ENB) (Murase & loka 2013; Fraija 2014a; Taboada 2010; Murase 
et al. 2014; Waxman 2013; Razzaque 2013; Murase et al. 2013) and 
to explain the recent detections of TeV- PeV neutrinos by IceCube 
(IceCube Collaboration et al. 2013; Aartsen et al. 2014). 
Measurements of HE neutrino properties such as flavor content 
would be involved with new physics if a deviation of the standard 
flavor ratio were observed (Learned & Pakvasa 1995; Athar et al. 
2000; Kashti & Waxman 2005; Mena et al. 2014). The neutrino fla¬ 
vor ratio is expected to be at the source, : <(>2^=1 : 2 : 0 

and on Earth (due to neutrino oscillations between the source and 
Earth) 4>l^ : : 1 : 1 and 4>l^ : : 4>% = l : 1.8 : 

1.8 for neutrino energies lesser and greater than 300 TeV, respec¬ 
tively (Kashti & Waxman 2005). Also measurement of a non-zero 
013 mixing angle coming from astrophysical sources could be rel¬ 
evant to clarify the neutrino mass hierarchy as well as CP violation 
searches in neutrino oscillations (Nunokawa et al. 2008; Bandy- 
opadhyay & et al. 2009; Forero et al. 2012). 

As known, neutrino properties are modified when they propagate in 
a thermal and magnetized medium. A massless neutrino acquires an 
effective mass and an effective potential. The resonant conversion 
of active neutrino from one flavor to another {v^ 44 tz.,.) due to 

thermal and magnetized medium has been explored in many astro- 
physical contexts and has had relevant consequences in the dynam¬ 
ics of them (Wolfenstein 1978a; Ruffert & Janka 1999; Goodman 
et al. 1987; Volkas & Wong 2000; Dasgupta et al. 2008; Erdas & 
Isola 2000; D’Olivo & Nieves 1996b, 1994, 1996a; D’olivo et al. 
1992; D’Olivo et al. 2003; Nbtzold & Raffelt 1988; Enqvist et al. 
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1991). For instance, Fraija (2014b) has showed that the effect of 
magnetic field in the dynamics of the fireball evolution of GRB 
was to decrease the proton-to-neutron ratio aside from the number 
of multi-GeV neutrinos expected in a neutrino detector. 

Neutrino oscillations in vacuum and by matter effects in the failed 
GRB framework (along the jet and envelope of the star) have been 
examined by many authors (Mena et al. 2007 ; Razzaque & Smirnov 
2010; Sahu & Zhang 2010; Osorio Oliveros et al. 2013; Fraija 
2014a) and although these authors have studied the oscillations on 
the surface of the star due to its envelope, the effect of thermal and 
magnetic field plasma generated on internal shocks has not been 
explored. In this paper we calculate the effect of the magnetized 
and thermal shocked plasma on neutrino oscillations and then we 
estimate the flavor ratio on Earth. The organization of the paper 
is as follows: In section 2, we show a brief description of internal 
shocks. In section 3, we derive firstly, the neutrino effective poten¬ 
tial for mw ^ as a function of the magnetic field, tempera¬ 
ture, angle (between the neutrino propagation and magnetic field) 
and chemical potential and secondly, the neutrino effective poten¬ 
tial produced by the envelope of the star. In section 4 we derive the 
resonance condition, the flip probability for two and three-neutrino 
mixing and the flavor ratio expected on Earth, and in section V we 
discuss our results. We hereafter use Qx = Q/KF in c.g.s. units 
and k=ft=c=l in natural units. 


2 DESCRIPTION OF INTERNAL SHOCKS 


One of the most prosperous theory to explain the prompt emission 
and the afterglow in successful GRBs is the fireball model (Zhang 
& Meszaros 2004; Meszaros 2006). A GRB is considered success¬ 
ful when the jet drills inside the progenitor and breaks through the 
stellar envelope, otherwise it is taken into account as a failed GRB. 
When the jet encounters the stellar envelope two shocks are in¬ 
volved: an outgoing, or forward, shock (Rees & Meszaros 1994; 
Paczynski & Rhoads 1993) and another one that propagates back 
decelerating the ejecta, the reverse shock (Meszaros & Rees 1994; 
Rees & Meszaros 1994). The jet dynamics is mainly dominated by 
the jet head, which is controlled by the ram pressure balance be¬ 
tween the reverse and forward shock. If the luminosity (Lj) is low 
enough and/or the density of the stellar envelope is high enough, 
then the hydrodynamic jet is collimated and internal shocks might 
occur inside the progenitor (Mizuta & loka 2013; Bromberg et al. 
2011; Murase & loka 2013). In this model, inhomogeneities in 
the jet lead to internal shell collisions, higher shells (Th) catching 
slower shells (Tj). The kinetic energy of ejecta is partially dissi¬ 
pated via these internal shocks which take place at a distance of 
rj = 2r^ ti, < i?*, where is the variability time scale of the 
central object, F ~ is the bulk Lorentz factor of the propa¬ 

gating shock and i?* is the radius of the progenitor’s stellar sur¬ 
face. The constraint Vj < i7* gives rise to those shocks inside 
the star. The physical width of the internal shock is lower by a 
factor of r, i.e. Arj = r/„. These internal shocks are expected 
to be collisionless, so that particles may be accelerated. In inter¬ 
nal shocks the total energy density U = 1/(8 vr trip) F”^ Lj 
is equipartitioned to generate and/or amplify the magnetic field 
eb = Ub/U = {B^/8^)11 (Piran 2005) and to accelerate parti¬ 
cles Ce = U^/U, where mp is the proton mass. Then, the magnetic 
field generated at the shocks is written as 


B' = e 


1/2 -p — 2 ^ 1/2 


( 1 ) 


It is important to say that the strength of the magnetic field falls 
out of the shocked region achieving some Gauss and although its 
direction might be random, it is mostly transverse to the jet di¬ 
rection (Razzaque & Smirnov 2010). From the causality condi¬ 
tion, the coherence length of such magnetic field is only the or¬ 
der of As ~ tv On the other hand, electrons are accelerated up 
to ultra-relativistic energies and then are cooled down rapidly in 
the presence of the magnetic field, producing the prompt emis¬ 
sion by synchrotron radiation. The opacity to Thomson scattering 
is Lj and photons thermalize at a black body 

temperature with peak energy given by (Razzaque et al. 2004) 


rp! ^ rl/4p —Ij-—1/2 

± ^ ? 
' TT 


( 2 ) 


where ctt is the Thompson cross section. 

Protons are also accelerated and cooled down in internal 
shocks via electromagnetic (synchrotron radiation and in¬ 
verse Compton (IC) scattering) and hadronic (proton-photon 
and proton-proton interactions) channels. Proton-photon and 
proton-proton interactions take place when accelerated protons 
interact with thermal keV photons (eq. 2) and proton density 
at the shock, n/ = l/(8 7rmp) F“^ Lj (Meszaros & Wax- 
man 2001). In both interactions HE charged pions and kaons 
are produced; p + ■y/p X jK^, and subsequently 

neutrinos tt"*’ —>■ pF + —>■ + Ve + Vp + and 

7r“ —> p~ + i/fj, —^ e~ + Ve + 1^11 + Vp. In this approach, the 
neutrino created by these processes will lie in the TeV - PeV 
energy range (Murase & loka 2013; Fraija 2014a; Razzaque & 
Smirnov 2010). 


3 NEUTRINO EFFECTIVE POTENTIAL 

In this section we are going to compute the neutrino effective po¬ 
tential due to the magnetized and thermal shocked plasma, and the 
envelope of the star. 


3.1 Magnetized and thermal plasma 

Recently, Fraija (2014b) derived the neutrino self-energy and ef¬ 
fective potential up to order mj/? at strong, moderate and weak 
magnetic field approximation as a function of temperature, chem¬ 
ical potential and neutrino energy for moving neutrinos along the 
magnetic field. In this subsection, we will calculate the neutrino ef¬ 
fective potential at the moderate and weak magnetic field limit for 
any direction of neutrino propagation. Therefore, following Fraija 
(2014b) we will show the equations that are more relevant for de¬ 
ducing the neutrino effective potential. 

The neutrino effective potential is calculated by means of the dis¬ 
persion relation 

14// = fco - |k| , (3) 

where k is estimated through the neutrino field equation in a 
medium (Nbtzold & Raffelt 1988; Enqvist et al. 1991) 

[k - E(fc)]t/)L = 0 , (4) 

and ai{k) = TZ{a\\k\\ -F apk±_ -F bit -F cb)C is obtained from the 
real part of its self-energy diagram. Here feJi* = (fc°, k^) and = 
(fc^, k^) are the momentum along and perpendicular to the mag¬ 
netic field, respectively, stands for the 4-velocity of the center- 
of-mass of the medium given by = (1, 0), 7?. = |(1 + 75 ) 
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and C = 1(1 — 75 ) are the projection operators and a, b, and 
c are the Lorentz scalars which are functions of neutrino energy, 
momentum and magnetic field. These scalars are calculated from 
the neutrino self-energy due to CC and NC interactions of neutrino 
with the background particles. The effect of the magnetic field is 
introduced through the 4-vector which is given by = (0, b) 
(Fraija 2014b). Using the Dirac algebra and from the dispersion 
relation (eq. 3), the neutrino effective potential can be written as 


where ?7 f(p • u) contains the distribution functions of the particles 
in the medium which are given by: 

, , 0{p ■ u) d(—p ■ u) 

tjfIp ■ u) - j;'. (12) 

where /3 and pe are the inverse of the medium temperature and 
the chemical potential of the charged lepton. By evaluating eq. (6) 
explicitly we obtain 


14// = b — c cos ip — a_L|k| sin^ p, 


(5) 


where p is the angle between the neutrino momentum and the mag¬ 
netic field vector. Otherwise, the effective potential that is applica¬ 
ble to the neutrino oscillations in matter is 14// = 14 — 14 ,t which 
depends only on electron density (Wolfenstein 1978b; D’olivo et al. 
1992). For that reason, although the one-loop neutrino self-energy 
comes from three parts; the lU-exchange, .Z-exchange and tadpole 
(Babaev 2004; Erdas et al. 1998; Sahu et al. 2009a,b), we will only 
consider the neutrino effective potential due to charged currents 
E(fc) = Ew(fc). We will use the finite temperature field theory 
formalism and the Schwinger’s propertime method to include the 
magnetic field (Schwinger 1951). From the W-exchange diagram 
(see fig. 1), the self-energy can be explicitly written as 


-iE(fc) = 7^ 


/ 


d*p 




( 6 ) 


where the weak coupling constant, is 

the W-boson propagator that in unitary gauge can be written as (Er¬ 
das et al. 1998; Sahu et al. 2009b) 





, 3 ie 


(7) 


here mw is the W-boson mass, Gf is the Fermi coupling constant, 
g^'' is the metric tensor and is the electromagnetic field tensor. 
From eq. (6), SKp) is the charged lepton propagator which is split in 
two propagators; one in presence of an uniform background mag¬ 
netic field (Slip)) and the other in a magnetized medium (Sj (p)), 
then it can be written as 


Se{p) = Slip) Sj(p). 


( 8 ) 


We can express the charged lepton propagator in presence of an 
uniform background magnetic field as 

roo 

iS°(p)=l s)ds, (9) 

JO 

where the functions ^(p, s) and G{p^s) are written as 

$(p, s) = is{pl - ml) - is[pl -t- , 

G(p, s) = sec^ 2[A+ *^75 

-|-m/(cos^ 2 — io-f sin 2 cos 2 :)], (10) 

where m; is the mass of the charged lepton, p^ = Po — P 3 > 
p3 = Pi + pi are the projections of the momentum on the mag¬ 
netic field direction and 2 = eBs, being e the magnitude of 
the electron charge. Additionally, the covariant vectors are given 
as follows, — Pfi — sin^ z{p ■ u — p ■ b bf , = 

sin 2 cos 2 (p ■ u bfj_ — p - b ufj, and erf = -)^bh . The other term in 
eq. (8) (due to magnetized medium) is given by (D’Olivo & Nieves 
1996a) 

/ OO 

e*(^’'*)G(p, s) ds , (11) 


ReT,{k) = TZ [aA_kF-\-bii +d)\C, (13) 

where the Lorentz scalars are given by (Fraija 2014b) 


ax = - - 


V2Gi. 


eB 

27r^ 


Ev^irie - rie) + - nf) 

H 


roo 

/ (2 
Jo „_n 


^n,o) 


(/e,n “1“ fe,n) 


(14) 


b = V2 Gf 

eB 


1 H- f - I (rie - ne) + (nf - nf) 

TTLt, 




n-2rn?_ 


rco f f 

■ J dp3 'Y2 (2 — Snfi)Eu^ I -I- 


2P] 1 I 3" /e.n) 


(15) 


and 

c = \J2Gf 
eB 


q 






1-^ (nf-nf)-^(ne-ne) 


j dp3 'Y2 (2 “ ^n,o)Ev^ I ‘^",0 


( „ 3 ml 


• (16) 

Here the electron number density and electron distribution function 

are 

R OO 

ne(M,T,B) = |^y](2-5„.o) 

n=0 

r dp3 

(17) 

Jo -p 1’ 


and 

fiEe,n,p) = + 1 ’ 


(18) 


respectively, with fe,n{p,T) = fe,n{—p,T) and Ee,n = 
yM + ml H with H = 2neB. Solving the integral terms in 
eqs. (14), (15) and (16) and replacing them in eq. (5) we calculate 
the neutrino effective potential for two cases: the moderate and the 
weak magnetic field limit. 


3.1.1 Moderate Magnetic field limit 


In the moderate field approximation {B j Be ^ 1), the Landau lev¬ 
els are discrete and can be described by sums (X)n 2, 3 

..). In this regime, the neutrino effective potential is written as 




Gf mjB 
m? Eu 


-4 


mfx me 


^(—1)‘ sinha; [Em — Gm cosp] 

i=0 

OO 

^( —1)* COShQ; [Jm — Hm COS <p] 


(19) 
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where ai = I3^{1 + 1 ) and the functions F^, Gm, Jm, Hm are 
written in the appendix A. It is worth noting that as the magnetic 
field decreases the effective potential will depend less on the Lan¬ 
dau levels. 


3.1.2 Weak Magnetic field limit 


In the weak field approximation {B j <C 1), all levels are full and 
overlap each other. In this regimen, sums over the Landau levels can 
be described and approximated by an integral —>■ f dn, then 
the effective potential does not depend on the Landau levels. The 
potential in this regimen can be written as 




eff,is{w) 


\/2 Gf tti^B 

TT^ Be 


^(- 1 )' sinha; [F^, - cos^s] 


2 j-i oo 

— 4 :^^ -^ y^(—1)* cosha; [Ju, — Jfu) cosipl (20) 


mf.r me 
^ 1=0 


where the functions F^,, Gw, Jm, are shown in the appendix A. 


3.2 Density profiles of envelopes 


is the mass difference (Giunti & Chung 2007), U is the two- 
and three-neutrino mixing matrix (see appendix B, eq. B9), Veff,k 
is the neutrino effective potentials calculated in section 3 (for fc=is 
and ss) and Ei, is the neutrino energy. We hereafter use the first 
and second line for two- and three-neutrino mixing, respectively, 
as written in eq. (25). From the conversion probabilities, we obtain 
that the oscillation lengths are 


Lres=4nEA cos29)2 + (5,m2 .in 29)2 - 

with the resonance conditions 

5m^ cos 29, 

Sm^2 cos 26^13. 


2 X 10 EvVeff.k 


(27) 


In addition to the resonance condition, the dynamics of the tran¬ 
sition from one flavor to another must be determined by adiabatic 
conversion through the adiabaticity parameter (Mohapatra & Pal 
2004) 


7 = -- 

2E„ 


1 


1 dVeff.k I 
V.tf.k dr \r 


Sm? sin 26 tan 26, 
5m|2sin26i3 tan 2613 , 


(28) 


Models of density distributions in CCSNe have been widely 
explored (Bethe & Pizzochero 1990; Chevalier & Soker 1989; 
Woosley et al. 1993; Shigeyama & Nomoto 1990). We will use two 
models with density profiles p oc and p oc Explicitly, 

the first model corresponds to a polytropic hydrogen envelope 

pi(r) = 4.0 X 10~® gcm“^, (21) 

and the second model is a power-law fit with an effective polytropic 
index rie// = 17/7 as done for SN 1987A (Chevalier & Soker 
1989) 

p2(f) = 3.4 X 10~^ g cm“® 

f ; 10^° ® cm < r < n = 10^^ cm 

\ {Ri,/(r - R*Y/{rb - R^^ ; r > n . 

( 22 ) 

In both cases, from the number density of electrons Ne = 
Na p(r) Ye, the neutrino effective potential can be written as 

Veff.as = V 2 GpNe , (23) 

where Na = 6.022 x 10^®p“^ is the Avogadro’s number. Ye =0.5 
is the number of electrons associated per nucleon and p(r) is given 
by eqs. (21) and (22). 


4 NEUTRINO RESONANT OSCILLATIONS 

When neutrino oscillations take place in matter, a resonance could 
occur that would dramatically enhance the flavor mixing and could 
lead to a maximal conversion from one neutrino flavor to another. 
This resonance depends on the effective potential and neutrino os¬ 
cillation parameters. The equation that determines the neutrino evo¬ 
lution in matter in the two and three-flavor framework is (Fraija 
et al. 2014) 

U ■ + diag{Veff,k,0) , (24) 

AiYji/ 

where 

_ J {—5m?,Q) for two flavors, 

\ (—(Jmli, 0 , Sm 32 ) for three flavors, 


with 7 3 > 1 or the flip probability given by 

= (29) 

By considering that the flux ratio of ~ 2Ni,,, ~ 2Nj2^ 

is created in the internal shocks, neutrinos firstly oscillate in matter 
due to the magnetized and thermal plasma and secondly oscillate 
to the star envelope. In vacuum, after neutrinos have left the star, 
they start oscillating to the Earth. Hence, from these three effects: 
internal shocks, envelope of the star and vacuum, the flavor ratio 
expected on Earth will be 



where the probabilities P*j are derived in appendix B. 

The best fit values of the two neutrino mixing are: Solar 
Neutrinos: 5m^ = (S.fiti ®) x 10“®eV^ and tan^ 6 = 

0.427lQ'Q29(Aharmim & et al. 2011), Atmospheric Neutrinos: 
5m^ = (2.11q' 4) X 10“® eV^ and sin^ 26 = I.OIq qy (Abe 
& et al. 2011) and Accelerator Neutrinos: 5m^ = 0.5 eV^ and 
sin^ 6 = 0.0049 (Zeitnitz 1994; Athanassopoulos & et al. 1996, 
1998) . Combining solar, atmospheric, reactor and accelerator pa¬ 
rameters, the hest fit valnes of the three neutrino mixing are, 
for sin^g < 0.053 ; = ( 7 . 4 l/(Q' 4 g) X 10~® eV^ and tan^ 612 = 

0.446to.S29 and, forsinfg < 0.04 : = (2.lto;|) x 10-3 eV^ 

and sin^ 023 = 0.50(tg Qg 3 , (Aharmim & et al. 2011; Wendell & 
et al. 2010). 


5 RESULTS AND CONCLUSIONS 

In this analysis we have considered HE neutrinos created in the 
energy range of 100 GeV^ E^ <100 TeV (Murase & loka 
2013; Eraija 2014a; Razzaque & Smirnov 2010) and also we 
have assumed (in the CCSNe-GRB connection) progenitors such 
as Wolf-Rayet (WR) and blue supergiant (BSC) stars with radii 
Ri, = 10 ^^ cm and R^, — 3 x 10 ^^ cm, respectively, with 
formation of jets leading to internal shocks inside of them. 

In internal shocks, energy is equipartitioned to generate and/or 
amplify the magnetic field and to accelerate particles. Electrons 
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and protons are expected to be accelerated in these shocks, and 
after to he cooled down hy synchrotron radiation, inverse Compton 
and hadronic processes (p 7 and p-hadron interactions). Photons 
produced hy electron synchrotron radiation are thermalized at 
keV energies and serve as targets for production of HE neutrinos 
through K^, and decay products in the proton -7 
proton-hadrons interactions. Therefore, this plasma is endowed 
with a magnetic field and made of protons, mesons, electrons, 
positrons, photons and neutrinos. 

First of all, we consider those internal shocks that take place 
inside progenitors (tj < R^,), as plotted in fig. 2. In this figure, 
we show the contour lines of hulk Lorentz factors and variability 
time scales for different internal shock radii. In these plots we 
observe that for a typical value of variahility in the range of 
10“^ s sC ^ 1 s, the values of hulk Lorentz factors are F ^ 34 
for a WR (left-hand figure) and F < 163 for a BSG (right-hand 
figure). Taking into account internal shocks at rj = 10^°'® cm 
(left-hand figure) and rj — cm (right-hand figure) we can 

see that the physical width of the internal shocks is restricted to 
Avj ^ 1.8 X 10^° cm and Avj ^ 2.5 x 10^^ cm, respectively. 
Once obtained the values of F and tj, for internal shocks inside 
the progenitor we compute the range of values associated to the 
magnetic field and temperature of the plasma, as shown in figs. 3 
and 4, respectively. We plot the contour lines of the magnetic fields 
(fig. 3) and thermalized photons (fig. 4) for values of luminosity 
in the range 10^® erg/s ^ ^ 10®^ erg/s. Colors in light- and 

dark-gray backgrounds represent the regions of a WR and a BSG, 
respectively. 

In figs. 3 and 4 one can see that the values of magnetic field and 
thermalized photons lie in the ranges 10 ® G s; B 10“ G and 
0.1 keV ^ T ^ 30 keV, respectively. It is important to clarify that 
the magnetic field amplified in the internal shocks falls out of them 
to the magnetic field endowed by the progenitor (black hole (BH) 
or magnetar) (Razzaque & Smirnov 2010). 

Following Fraija (2014b) and taking into account that the range 
of neutrino energy considered is larger than the W-boson mass 
ml ^ m^), we have obtained the neutrino effective potential 
up to an order in the moderate below B/Bc ~ 10“® and 
weak BjBc ~ 10 “^® regime as a function of the observable 
quantities in the internal shocks: thermalized photons, magnetic 
field, neutrino energy and angle between the direction of neutrino 
propagation and the magnetic field. We plot the neutrino effective 
potential in both limits (moderate and weak field limits) as 
shown in fig. 5. The neutrino effective potential at moderate limit 
(figures above) and weak limit (figures below) are plotted for a 
magnetic field in the range of 10 “® Be < B < 10 Be and 
10 "®® Be < B < 10 Be, respectively. In both cases, we use 
the values of temperature T= (20, 24, 27 and 30) keV, angle ifi= 
(0°, 30°, 60° and 90°) and the neutrino energy E^=10 TeV. The 
neutrino effective potential at the weak limit is smaller than at 
the moderate limit. It is worth mentioning that in the range of the 
magnetic field considered, the contribution of Landau levels to 
the effective potential at the moderate-field limit is not significant 
due to A„ Ki{aiX„) ~ 0 for A„ = ^/T+~2dlB~/~^ and 

ai = Pme(l + 1). From these plots one can observe that the neu¬ 
trino effective potential is positive, therefore due to its positivity 
(Vef f,is(k) > 0 ) for k= m and w, neutrinos can oscillate resonantly. 
From the resonance condition (eq. 27) and the neutrino effective 
potential at moderate (eq. 19) and weak (eq. 20) limit, we plot the 
contour lines of temperature and chemical potential as a function 
of neutrino energy for which the resonance condition is satisfied, 
as shown in figs. 6 and 7, respectively. From these figures, one can 


Energy 

(TeV) 

Solar 

^res 

Atmosph. 

(cm) 

Accelerat. 

Three flav. 

10“® 

4.8 Xl0®° 

1.2 xlO® 

7.1 xlO'^ 

2.6 XlO® 

0.5 

2.4 xlO®® 

5.9 xlQi® 

3.6 XlO® 

1.3 XlO®® 

10 

4.8 xlO®® 

1.2 xlO®® 

7.1 XlO®® 

2.6 XlO®® 

10® 

4.8 xlO®"® 

1.2 xlO®® 

7.1 XlO®® 

2.6 XlO®® 


Table 1. Resonance lengths of HE neutrinos for the best fit parame¬ 
ters of the two- and three-neutrino mixing. 

see that temperature is a decreasing function of chemical potential 
and neutrino energy. As neutrino energy increases, temperature 
decreases steadily. Considering the values of neutrino energy (B^ 
= 100 GeV, 500 GeV, 10 TeV and 100 TeV) and (fi = 90°, we 
see that the temperature and chemical potential are in the range 
10 keV to ~100 keV and 60 eV to 50 keV, respectively. For 
instance, taking into account a neutrino energy of 10 TeV, from 
fig. 6 we can see that temperature lies in the range 22.2 to 15.3 
keV for solar, 23.4 to 16.2 keV for atmospheric, 40 to 26.4 keV 
for accelerator and 28.3 to 15.4 keV for three-neutrino parameters, 
and as shown in fig. 7, temperature lies in the range 18.3 to 14.1 
keV for solar, 19.8 to 15.1 keV for atmospheric, 32.1 to 21.9 keV 
for accelerator and 24.5 to 18.3 keV for three-neutrino parameters. 
In addition, we have obtained the resonance lengths which are 
shown in table 1. As shown in this table, the resonance lengths 
lie in the range Les ~ ( 10 ®® to 10 ®®) cm, hence depending on 
the progenitor associated and the oscillation parameters, neutrinos 
would leave the internal shock region in different flavors of 1 : 2 : 0 . 
For instance, taking into account the parameters of three-neutrino 
mixing, neutrinos with energy less than 0.5 (10) TeV will oscillate 
resonantly with a resonance length equal or less than the radius 
of the progenitor, either a WR or BSG. Considering parameters 
of accelerator experiments, neutrino energy around 100 TeV will 
oscillate in a BSG star before leaving it. 

As the dynamics of resonant transitions is not only deter¬ 
mined by the resonance condition, but also by adiabatic con¬ 
version, we analyze the flip probability (eq. 29) to find the 
regions for which neutrinos can oscillate resonantly. First of 
all we derive the neutrino effective potential as the function 
of magnetic field dVeff/dr = dVeff/dB x dB/dr, and 
assume that at internal shocks (10®® ® cm for WR and 10®® 
cm for BSGs), magnetic fields change a 10% of any variation 
around the radius shock scale. For instance, for a WR star, 
dBjdr = 0.1 X 10"'®Bc/10®®'®cm = 6.99 x 10“® Gauss/cm. 
We plot the flip probability as a function of neutrino energy for two 
and three flavors (fig 8 ). We divide each plot of flip probability in 
three regions in order to analyze the whole range of probabilities: 
less than 0.2 (P^. ^ 0 . 2 , a pure adiabatic conversion), between 
0.2 and 0.8 ( 0.2 < P-, < 0.8 represents the transition region) 
and greater than 0.8 (P.y 0.8 is a strong violation of adiabatic- 
ity)(Dighe & Smirnov 2000). In fig. 8 , we use two flavors: solar 
(left-hand figure above), atmospheric (right-hand figure above), 
accelerator (left-hand figure below) and three flavor (right-hand 
figure below). When we use solar parameters, a pure adiabatic 
conversion occurs in a WR (BSG) star for neutrino energies of 
less than 10®® ® (10®®“) eV and 10®®'® (10®®“®) eV which are 
endowed with B = 10“'®Be and B = 10“®®Bc, respectively. 
Considering atmospheric parameters, only a pure adiabatic conver¬ 
sion takes place in a WR (BSG) star for neutrino energies less than 
10®®.6 ( 1014 . 8 ) gy ;^q 14.5 JQIS) gy ^^.g gndowed 

with B = 10“'®Be and B = 10“®®Be, respectively. Taking 
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into account accelerator parameters, a pure adiabatic conversion 
happens in a WR (BSG) star for neutrino energies of less than 
IQii-S (iQis i) eV and > 10^®) eV with B = 

and B = respectively and once again considering three 

neutrino mixing, a pure adiabatic conversion occurs in a WR 
(BSG) star for neutrino energy of less than (10^^ ®) eV 

and (10®-®) eV with B = lO'^Bc and B = 

respectively. Higher energies to those considered are found in 
regions of transition and/or those prohibited. 

In addition, we have studied the HE neutrino oscillations from 
the neutrino effective potential generated in the star envelope (eq. 
23), as shown in fig. 9. From the resonance condition (eq. 27), we 
obtain the contour plots of radius as a function of neutrino energy. 
One can see that for neutrino energy in the range 100 GeV< < 
100 TeV the radius lies in the range 10^°'® cm < r < 10®^ ® 
cm. The flip probability for neutrino oscillations in the envelope 
of a star was studied by Fraija (2014a). The author has plotted this 
probability as a function of neutrino energy for density profiles 
[A] (eq. 21) and [B] (eq. 22) and neutrino oscillation parameters. 
From this analysis, Fraija (2014a) showed that neutrinos can 
oscillate depending on their energy and the parameters of neutrino 
experiments, obtaining that neutrino with energies above dozens 
of TeV can hardly oscillate. 

Finally, considering a flux ratio ~ ~ ~ 2Nji^, we 

estimate the neutrino flavor ratio coming from the surface of a WR 
and BSG to Earth, as shown in fig. 10. In this estimation, we take 
into account the contribution of thermal and magnetized plasma 
at moderate and weak B limit generated by internal shocks; at 
10 cm (second panel), 10 cm (upper panel), 10 cm (bottom 
panel) and 10 ^®'® cm (third panel), the effective potential due to 
the envelope of star and oscillation neutrinos in vacuum, due to 
the path up to Earth. In this figure we take into account two values 
of 013 mixing angle, 2 ° (left column) and 11 ° (right column). 
As shown, one can observe that a nonsignificant deviation of the 
standard ratio ; 1 : 1 : 1 ) is expected, 

less than 10 % for @13 = 11° and only 2 % for 613 = 2°. In 
addition, we plot the neutrino flavor ratio expected on Earth as a 
function of neutrino energy when the magnetic field is oriented 
to different angles 0° 75° concerning neutrino direction, 

as shown in figs. 11 and 12. In fig. 11 we consider the neutrino 
effective potential at the moderate-field limit and internal shocks at 
rj = 10 ®® cm with a physical width Ar^ — 2 x 10 ®® cm and in fig. 
12 , we consider the neutrino effective potential at the weak-field 
limit and internal shocks at rj = 10 ®*®'® cm with a physical width 
Arj = 1.5 X 10®® cm. From both plots, we can see that although 
the neutrino flavor ratio changes at different angles, distances of 
internal shocks, strength of magnetic field (moderate and weak 
limit) and neutrino energy in the range 10®® eV ^ ^ 10®'* eV, 

this flavor ratio expected on Earth lies between 0.98 and 1.02, 
hence we can conclude that the directionality of magnetic fields 
does not affect our results. Although currently neutrino oscillations 
can hardly be detected, new techniques in the near future will allow 
us to perceive these oscillations and put limits on the neutrino 
mixing angles. Finally, it is worth noting that the estimated values 
of the bulk Lorentz factor, in particular those relying on variability 
time measurements, are only raw approximations, and variations 
by a factor of a few cannot be ruled out by existing data. 
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APPENDIX A: EFFECTIVE POTENTIAL 


The functions of the neutrino effective potential at moderate mag¬ 
netic field limit are 

Frr. = \l + 2^-\Kdol) + 2ff\Jl + ^-\Kdm\f 
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and at weak magnetic field limit are 
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where Xf = l-|-2r!,B/ Be, Ki is the modified Bessel function of 
integral order i, q; = I3^,{1 + 1) and ct; = pmeil + 1). 


APPENDIX B: PROBABILITIES 


The flavor ratio at the internal shocks and the envelope of the star 
are 



(Bl) 


(B2) 


respectively Here ss and is are the envelope of star and internal 
shocks, respectively. The different neutrino probabilities can be 
written as (Gonzalez-Garcia & Nir 2003; Gonzalez-Garcia & Mal- 
toni 2008) 
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where 
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Figure 1. W-exchange diagram of one-loop contribution to the neutrino self-energy in a magnetized medium. The dashed line represents the electron propagator 
e~ {p), the solid line corresponds to the electron neutrino propagator Ue{k) and the wiggly line is the W-boson propagator W'^{q). 




Figure 2. Contour lines of variability time scale and bulk Lorentz factor (P) as a function of the distance of the internal shocks (rj) for which these shocks 
take place inside the progenitors. We have considered progenitors such as WR (left-hand figure) and BSG (right-hand figure) stars. 


sin 2013 ,m = 
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and 



Here are given by 
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2007; Bahcall 1989) 

= \< Vp{t)\Va{t ^ Q) >\ 
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3>i \ J 

where the neutrino mixing matrix Uij is given by Gonzalez-Garcia 
& Nir (2003); Akhmedov et al. (2004); Gonzalez-Garcia (2011) 
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where 
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(B7) 


where the neutrino effective potentials Veff,k are given in section 3. 
In vacuum, the flavor ratio (between the surface of the star and the 
Earth) is affected by the oscillation probabilities (Giunti & Chung 


( C13C12 S12C13 S13 \ 

—S12C23 ~ S23S13C12 C23C12 — S23S13S12 S23C13 1 , (B 9 ) 

S23S12 — S13C23C12 —S23C12 — S13S12C23 C23C13. J 

Here Sij = sin0ij and dj = cos 6ij and we have taken the Dirac 
phase (5 = 0. Taking into account the effect of internal shocks, 
envelope of star and the vacuum, the probabilities are given by 

Pi*i = 0.82 Pii+0.55 P 21 -h 0.19 P 31 

Pi *2 = 0.82 P 12 + 0.55 P 22 + 0.19 P 32 

Pi *3 = 0.82 Pi 3 + 0.55 P 23 + 0.19 P 33 

P21 = - 0.51 Pii + 0.51 P21 + 0.69 P31 

P 2*2 = -0.51 P 12 + 0.51 P 22 + 0.69 P 32 

P 23 = -0.51 P 13 +0.51 P 23 +0.69 P 33 

P3*i = O.28P11-I— 0.66 P31 + 0.69 P31 

P 32 — 0.28 P 12 H—O. 66 P 22 + 0.69 P 32 

P 33 ~ O.28P13-I— 0.66 P23 + 0.69 P33 (BIO) 
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Figure 3. Contour lines of variability time scale (ti,) and bulk Lorentz factor (F) as a function of magnetic field {B') generated at internal shocks (eq. 1) inside 
a WR (light gray background) and a BSG (dark gray background). In the left-hand figure above (below) a luminosity of 10'^® (10^^) erg/s was used, whereas 
in the right-hand figure above (below) it corresponds to a luminosity of 10*^® (10^^) erg/s. 
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Figure 4. Contour lines of variability time scale and bulk Lorentz factor (F) as a function of synchrotron photons created in internal shocks and thermalized 
to a black body temperature (Tp (eq. 2). Once again we consider a WR (light gray background) and a BSG (dark gray background). In the left-hand figure 
above (below) a luminosity of 10“^^ (10^®) erg/s was used whereas in the right-hand hgure above (below) it corresponds to a luminosity of 10“^® (10^^) erg/s 
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Figure 5. Neutrino effective potentials at moderate (top) and weak (bottom) limits are plotted as a function of magnetic field for temperatures at keV energies 
(left) and different angles between the direction of neutrino propagation and magnetic field (right). 
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Figure 6. Contour lines of temperature and chemical potential as a function of neutrino energy for which the resonance condition is satisfied. We have used 
the neutrino effective potential at the moderate limit (eq. 19) and the best-fit values of the two-neutrino mixing (solar, top left; atmospheric, top right; and 
accelerator, bottom left) and three-neutrino mixing (bottom right). 
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Figure 7. Contour lines of temperature and chemical potential as a function of neutrino energy for which the resonance condition is satisfied. We have used the 
neutrino effective potential at the weak limit (eq. 20) and the best-fit values of the two-neutrino mixing (solar, top left; atmospheric, top right; and accelerator, 
bottom left) and three-neutrino mixing (bottom right). 
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Figure 8. The flip probability is plotted as a function of neutrino energy for a strength of magnetic field in the moderate (B=10~^ Be) and weak (B=l G) 
regime and at a distance of 10^® ® cm and 10^^ cm. We have used the best-fit values of the two-neutrino mixing (solar, top left; atmospheric, top right; and 
accelerator, bottom left) and three-neutrino mixing (bottom right). 



Figure 9. Contour lines of distance and neutrino energy as a function of neutrino oscillation parameters for which the resonance condition is satisfied. We have 
used the neutrino effective potential generated by the envelope of star (eq. 23) and the best parameters of neutrino oscillation for solar, atmospheric, accelerator 
and three flavors. 
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Figure 10. Neutrino flavor ratio expected on Earth as a function of neutrino energy when these are created on the surface of a WR (at 10^^ and cm 

for the upper and second panels, respectively) and a BSG (at and 10^^ cm for the third and bottom panels, respectively). We have used the neutrino 

effective potential at the moderate (second and bottom panels) and weak (upper and third panels) field limit for 0 i 3 =l 1 ° (left) and ^13 = 2 °. (right) 
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Figure 11. Neutrino flavor ratio expected on Earth as a function of neutrino energy when these are created on the surface of a BSG. We have considered that 
internal shocks take place at rj = 10^^ cm with a physical width Ar^ = 2 x 10^^ cm and the magnetic field is oriented to different angles 0° ^ ^ 75° 

concerning to neutrino direction. We have used the neutrino effective potential at the moderate-field limit (eq. 19) and 613 = 2°. 
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Figure 12. Neutrino flavor ratio expected on Earth as a function of neutrino energy when these are created on the surface of a BSG. We have considered that 
internal shocks take place at Vj = 10^^ cm with a physical width Ar^ = 2 x 10^^ cm and the magnetic field is oriented to different angles 0° ^ ^ 75° 

concerning neutrino direction. We have used the neutrino effective potential at the moderate-field limit (eq. 20) and ^13 = 2°. 
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